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Thermodynamics of classical frustrated spin chain at the ferromagnet-helimagnet 

transition point 
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(Dated:) 

Low-temperature thermodynamics of the classical frustrated ferromagnetic spin chain is studied. 
Using transfer-matrix method we found the behavior of the correlation function and zero-field sus- 
ceptibility at the ferromagnetic- helical transition point. It is shown that the critical exponent for 
the susceptibility is changed from 2 to 4/3 at the transition point. 



Lately, there has been considerable interest in low-dimensional spin models that exhibit frustration. One of them 
^ ] is the spin chain with the ferromagnetic interaction Ji of nearest neighbor (NN) spins and the antifcrromagnetic 
'■^ ■ next-nearest-neighbor (NNN) interaction J2, so called the ID F-AF model. Its Hamiltonian has a form 

^ : 

(N . H = .hJ2 s«-s„+i + ./2 J2 s«-s„+2 (1) 

(— I , where Ji < and J2 > 0. 

Q ■ This model is characterized by a frustration parameter a = -h/lJil- The ground state properties of the quantum 

D I 5 = 1/2 F-AF chain have been intensively studied last years [l|-|j]- It is known that the ground state of the model is 
H ■ ferromagnetic for a < 1/4. At a = 1/4 the ground state phase transition to the incommensurate singlet phase with 
I [ helical spin correlations takes place. Remarkably, this transition point does not depend on a spin value, including the 
C^ . classical limit s = 00. 

"^ ' Interesting question is the influence of the frustration on the low-temperature thermodynamics of the model es- 

I ; , pecially near the transition point a = 1/4. We study this problem for the classical version of model ([T]). At zero 
Ci ' temperature the classical model has long range-order (LRO) for all values of a: the ferromagnetic LRO at a < 1/4 and 
H ', the helical one at a > 1/4. At finite temperature the LRO is destroyed by thermal fluctuations and thermodynamic 
I , quantities have a singular behavior at T ^ 0. In particular, the zero-field susceptibility x diverges. For the ID 

^ Heiscnbcrg fcrromagnet (a = 0) x = 2 \Ji\ /3T^ [|. At < a < 1/4 the susceptibihty is x = 2(1 - 4a) \Ji\/3T^. 
This behavior of x is similar to that for the quantum s — 1/2 F-AF model [6|. The prefactor in x vanishes at the 
fj ' transition point indicating the change of the critical exponent. We focus our attention on the behavior of x at the 
transition point. 

The partition function Z of model ([T|) at a = 1/4 is 



Y[ / dfln exp <^ - ^{Sn ■ Sn+1 ~ -Sn ' Sn+2) > 



(2) 



where S'„ is unit vector, dfln is the volume element of the solid angle for ?i-th site, we put | Ji| = 1 and the periodic 
boundary conditions are proposed. 

Our further calculations are based on the transfer matrix method and we use a version of this method adapted to 
the model with NNN interactions by Harada and Mikeska in [3j • 

Following Ref.[7j we represent Z in a form 

N 



X 

C^ n=l 

where 



Z=ll I dnr.Ki0n-l,en;ipn) (3) 



/cos6^„_i +cos6'„ (cosy„_iCOsy„ +smy„_ismf^„cos.^„) \ 
K{9n-i,0n;(pn) = expl — — I (4) 

where 9n is the angle between Sn and Sn+i and tpn is the angle between components of Sn-i and Sn+i projected 
onto (XmYn) plane of the 71-th local coordinate system with the Z„ axis parallel to Sn- 
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Integrating Eq.Q over ifn we obtain Z in a form 

Z = \\l d9n sin 9nA{9n-i,9n) 

n Jo 

where 



(5) 



Ij, . /cos6'„_i +cos6'„ cos6l„_icos6i„ 3 , 
A(6'„_i,6'„) == -/o(-z)exp ( — — — ) (6) 



and /o(— z) is the modified Bessel function of 

sin0„_i sin^n 



z 



AT 
Let us consider an integral equation 



(7) 



Ai9i,92)M02) sin (^2) d92 = XcM^i) (8) 



where ^a{9) satisfy normahzation condition 

/ i>c.{e)^fi{e)sin9A9^6^,0 (9) 

Jo 

Eigenfunctions V'a and eigenvalues A^ can be chosen as real, since the kernel A{9i^ ^2) is real and symmetric. Then, 

^(^l,^2)=5I^"^"(^l)^"(^2) (10) 

a. 

Substituting Eq. (|1QP into Eq/([5]) we obtain in the thermodynamic hmit 

Z = \^ (11) 

where Aq is the largest eigenvalue of Eq. ([5]) . 

In the low-temperature limit the angles 9n are small and we can use the asymptotic expansion of the modified 
Bessel function 

/o(-z) = ^fl + ^ + 0(z-^)) (12) 

Then, we expand the expression in the exponent of the transfer matrix to the fourth order in 9i to obtain 



A(.„.2)^,/^ri+^Vxpf-(^^-^-)^-gM+(^^-^^r) (13) 

^ ' ' y 2tt9i02 V 26*16127 \ 8T 8T 96T ' ^ ^ 

We can neglect the term {9i — 02)^/96T as will be seen below. 
As a result integral equation ([8]) reduces to 



V il(^ + ^) ^^-^^-^Mo.)de2 = X.MO.) (14) 

The maximum of the expression in the exponent (saddle point) is at 62 — 9i (more exactly 62 — 9i — 91 + . . ., but 
it suffices to put 02 = ^i)- Near this saddle point we expand 4'a{92) as follows 

^-0(^2) = V^a(ei) + [92 - ei)C(^l) + ^^'~/'^ C(^l) + ■ ■ ■ (15) 

and 



l)^ 



, J2 / l7? — "1 C/o ^ (71 1(79 — (71 1 . . 

'/-a/i + ^-1 + ^-^^^7^ + --- (16) 



'1 



Let us introduce new scaled variables 



Now -0Q (9) -^ ipa (r) and 



92-91 = T^l'^x 
9i = T^/^r 



MO2) ^ i^air) + T'/'x^P'^{r) + ^^CM + 



(17) 



(18) 



rpi/6^ rpi/3^2 



(19) 



exp 



Summarizing all above we arrive at 



n-«^2)2 9l9l\ I X 

' — ^ exp — -— 



ST 



8T 






2 2^1/3^,4 



8 8 



(20) 



(^^^{r)+T^/^xi^'M + ^^€m) e 



-xVa-ri/^rVsZi^ 



AQV'Q(r) 

(21) 

At T — > 0, we can change the limits in the integral to [—00, 00], then only even powers in x gives contribution, so 
taking into account only terms up to T^'^ we obtain 



yl/3^4 J.1/3X yl/3^2 / , , 



8 2r2 y "" ' 2 \ ^ r 

After integration over x we obtain a linear differential equation 

yi/3^4 2.1/3 



/^j'2 



, Tdx 



— -^aWa 



2T 1 



+ ^) V-a + 4T4/3 (^V^ + iv^(0 - ^^a) = A„7^„ 



and, finally, 



with 



1 r* 



Cry - 



2r- A„ 
4T4/3 



(22) 



(23) 



(24) 



(25) 



Thus, we have got a Schrodinger equation for a particle with Z component of the angular momentum Zz = in 2D 
potential well U{r) — r'*/16. Normalization condition for ipaif) is 



/•CX) 

/ ipa{r)^f3ir)'2T:rdr ^ Sa,p 

JQ 



Numerical solution of Eq. (|24|) gives the following lowest eigenvalues (corresponding to the largest A): 

e„ = 0.9305; 3. 78; 7.44... 
As was shown in Ref.jTj the two-spin correlation function can be expressed by the following integral 



Si ■ S'l+r 



1 r / n 

— ^ / d9nsm9n n d9i sin 9iM&i)M&n) (01) B{9i)H{9i,9i+i)B{9n) . 
Aq Jo i^-^ V ^ 



where 



B{9)^ 



cos6'/2 sin 61/2 
-sin 61/2 cos6'/2 



(26) 



(27) 



(28) 



(29) 



H(eu9 



'l,f2 



m)(-^>"«^„%J«(« 



(30) 



and A{9i,92) is given by Eq.® with Io{—z) replaced by Ii{~i 
Using the asymptotic expansion of the Bessel function 



h{-z) = - 



we obtain 



where 






3T 



\ 9 ^ ' 9fl,l 



26/16/2 



^°(^-^^)-vS^^^p(-^ 



9i - 02)^ 0?0i 



8T 



(31) 



(32) 



(33) 



The matrix H{9i,d2) is not symmetric. Therefore, to calculate (Si ■ 5i+„) it is necessary to solve a pair of the 
integral equations 



/■TT 

/ H{9i,92)MS2)sm{92)d92 = V^M^i) 
Jo 

f H^{9i,92)Vc.{92)sm{92)d92 = r?a«a(0l) 

Jo 
where H^ {9i,92) is transposed matrix H{9i,92) and two-component vectors Ua and Va 



U2,a 



Va = 



^1, 
«2, 



satisfy orthonormality relations, 



ul{0)vp{9) sin (6*) de* = / C(^)",3(6') sin {9) d9 = (5„,/3 

/n 



Then, the matrix H{9i,92) can be represented as 

i^(^i,e2)=5Z^"^"(^i)^"(^2) 



(34) 

(35) 

(36) 

(37) 
(38) 



At smaU 9i, 92 Eqs.(l34]) and ([35]) reduce to 



^0(01,^2) 



1- 



3T 

2M^ 



Ml,a(6'2) + 6'iU2,a(6'2) sin(6'2)d6'2 = riaUi^a{9l) 



Ao{9i,92) -9iui^a{e2) + (l + ^^] U2A02) sm{92)d92 = r7a"2,a(^i) 



(39) 
(40) 



Integrating these equations near the saddle point similar to Eqs. (l2ip . (j22p . we get a pair of linear differential equation 



2T 
2T( 1- 



8 2r2 

2nl/3j,4 yl/3 



(41) 



2r2 



U2.a 



- 4T4/= 



'3 f „/' 



-"2,Q ~ l-l-"2,a j - 2r'*/^rMi,a = 77aU2,a 



'2,a+ ^W2,Q ^^2' 



(42) 



and, finally, 

" 1 ' +i 

where 



l,a - -'^l,a + — Wl,a + 77^1, a + ^^"2,0 = Ma^l.a (43) 

W2,Q - r'"2,a + T^"2,a " T^^l^a = ^J■aU2,a (44) 



'^^ = 4r4/3 (45) 

A few lowest eigenvalues of Eas. (|43|) and p4|) are 

/Ua = 1.4113; 1.83; 3. 98... (46) 

For Va similar procedure gives 



1,1 r r 

J, , j,^ 16 



"l,a "Vl^a + —'^l'" + Tg^I^" ~ o^2,a = A^Q^l.a (47) 



-«2,a - -«2.a + T7^2,a + ^rWl.a = Aia'y2,a (48) 

r 16 2 

It follows from Eas. p5)) - (PH) and (H7|) - (H51) that the functions v^ is connected with u^ by the relations vi^ — — Wi,q, 
V2.a = U2,Q and, therefore, normalization condition ([37]) transforms to 

/"OO 

r^^^ / {u2.aU2.j3 - ui^aUi.fj) rdr = (Sa_^ (49) 

Jo 

Using Eqs. (f38|) and (j37p we obtain the correlation function (|28|) in a form 

(5'i-5'i+„)=^2;rV^ (50) 

a 

where j/q = rja/Xa and 



At T^O 

and the correlation function becomes 

[Si ■ 5i+„) = ^ /^ exp[-2ri/3 (^, - eo) (n - 1)] (53) 

a 

According to Eg. ([53]) the correlation length ^ at T -> is 

^"2(Mo-£o)Ti/3-TV3 ^^^> 

Now we are ready to calculate the magnetic susceptibility at T — > 0, which is 



U= 7/^o(%2,a(0)sin(0)d0-T2/3 / Mr)u2Ar)rdr (51) 



— ^(5i.5i+„)_— (1 + 2^^-^)- — + ^^^^- 



x = ^1.(^-^iW-3t(^+21^t3^)-3t + 3^^;:::^ (^^) 



Now we see that /^ and (/Iq, — Eq) depends on the solutions of differential equations which are independent of T. 
So, the sum in x gives numerical constant 

y •'" = 3C (56) 

^ lla- So 



Therefore, the low-temperature susceptibihty behaves as 



C 



X- 



Ji 



1/3 



T4/3 



(57) 



Numerical calculations gives for the constant C the value C « 1.07. Thus, the critical exponent for the susceptibility 
at the transition point is 4/3 and that for the correlation length is 1/3. 
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